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The q state Potts model is one of the most intensively 
studied classical two-dimensional lattice models.!!© The 
model shows second-order phase transitionpjdjcn q < 4, 
and shows first-order transition when q > 4.aa The crit- 
ical phenomena of the case Jj <5 4 have been classified by 
the conformal field theory.t§El ) Since critical exponents 
are exactly known for q = 2,3, and 4, the Potts model 
has been used -as. a reference system for numerical fi- 
nite size scaling!! u (FSS) analyses. Recently Carlon and 
IgloitP showed tdiai. .the density matrix renormalization 
group (DMRG) ffiEjiLM jrfdeb has been applied to classi- 
cal lattice modelsJi3 , E£toE3 ) is applicable to the precise 
estimation of the critical exponents. They also suggest 
the applicability of DMRG to the first-order phase tran- 
sitiorJIlia) 

The phase transition of the two-dimensional q = 5 
Potts model is often referred as 'weak first-order tran- 
sition', because the latent heat L ^jO.0265 is especially 
small compared with the cases q > 6& To distinguish the 
weak first-order transition from the second-order tran- 
sition is a remaining problem in numerical analyses of 
statistical models, and therefore tow-dimensional. O-rj— 5 
Potts model has been numerically investigated Ji2Hi3'c3 
At present, it is possible to detect non second-order ten- 
dencies in the order parameter^ or the specific heat.Es' 
However, quantitative observation of the latent heat L 
has not been performed. In this short note we calculate 
the latent heat L of the q = 5 Potts model directly by 
subtracting the site energy of the ordered (= ferromag- 
netic) phase at the transition temperature from that of 
the disordered (= paramagnetic) phase. For this purpose 
we use the corner tmnafer matrix renormalization group 
(CTMRG) methodBBBEf 

The outline of CTMRG is as follows. Let us con- 
sider a square cluster of g-state Potts model whose lin- 
ear size is N. The cluster consists of four sub-clusters 
of the. linear size (N + l)/2, which are called as 'cor- 
ners'.® The partition function Zjq is calculated as the 



CTM_,BB ) and transform A {N+1)/2 into a renormalizcd 
one j4(at-|-i)/2, which is an m-dimensional matrix. Trace 

of pn = {A(n+i)/2) gives a precise lower-bound Zm 
for the exact partition function Zjy, where the differ- 
ence Z/v — Zn approaches to zero with increasing to. 
We can increase the linear size of the corner one by one 
using a recursive relation between A n and A n+1 , and ob- 
tain ^4(at+i)/2 for arbitrary N. The renormalized CTM 
A(n+i)/2 thus obtained successively gives the approxi- 
mate partition function Zjy. In addition to the partition 
function, we can calculate the site energy and the order 
parameter at the center of the square cluster. 

The partition function of the two-dimensional q = 5 
Potts model is given by 



rN 



(ij) 



(0.1) 



where N is the linear dimension of the finite size sys- 
tem, <7i is the five-state spin variable on the lattice point 
i, (ij) specifies the neighboring spin pairs, and 5((7i,<jj) 
is equal to unity only when cr^ = aj and zero other- 
wise. We consider both paramagnetic boundary con- 
ditions (PBC), where arbitrary spin configurations are 
allowed at the boundary, and ferromagnetic boundary 
conditions (FBC), where boundary spins are fixed to the 
same direction. In the limit of TV — » oo, the q = 5 Potts 
model shows first-order phase transition when the pa- 
rameter K is equal to 



K c = ln(V5 + 1) . 
At the transition point, the site energy 
U = 2{5(a i ,a j )) 



(0.2) 



(0.3) 



trace of pn — (A 



(JV+l)/2j 



where A 



(N+l)/2 



of the disordered (= high-temperature) phase U + differs 
from that of the ordered (= low-temperature) phase U~ 
by the latent heat L ~ 0.0256. The average of U + and 
U~ is simply expressed as Uq =p.J/2 + 1/2 v5, that is 

nil) 



is so called obtained by the duality relation! 



'corner transfer matdx' (CTM) which transfers column 
spins into row spinslf The point of CTMRG is that p^ 
can be regarded as a kind of density matrix in DMRG, 
because its trace is the the partition function Zjy. We 
can therefore apply the renormalization procedure to 
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In order to estimate numerically, we first calcu- 
late the site energy U^(m) and U^(m) at the center of 
square cluster of size N, respectively, by imposing PBC 
and FBC when K = K c . The number to is the matrix di- 
mension of the renormalized CTM A(jv+i)/2- It should 
be noted that U ± is equal to lim m _ >0O limAr^ oo [/^(TO,). 
Figure 1 shows U^(m) for m = 40,67 and 200 up to 
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Fig. 1. System size dependence of the site energy for both fixed 
and free boundary conditions. 



Fig. 2. The m dependence of U 00 (m). The 1/m extrapolation 
gives the latent heat L. 



N = 4000. We choose Uq as the zero of the site en- 
ergy. The m dependence of U^(m) is not negligible in 
the large ./V region, where the energy (m) approaches 
to U^(m) exponentially with respect to N. 

Figure 2 shows the m dependence of U^(m). There 
is m dependence because the energy cut-off introduced 
by the renormalization group transformation is compa- 
rable to the latent heat L. As shown in Fig. 2, U^(m) 
is nearly linear in 1/m, and therefore we may linearly 
extrapolate U^(m) to obtain = U^(oo). As a re- 
sult we obtain U+ = U - 0.014 and U~ = U + 0.013, 
where the obtained U + and U~ are slightly lower than 
the exact ones. We finally obtain the latent heat L = 
U~ — U + = 0.027, which quantitatively agrees with the 
exact one L ~ 0.0265. As far as we know, this is the 
first quantitative estimate of the latent heat by numeri- 
cal calculation. 
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